Charge accumulation at the boundaries of a graphene strip induced by a gate voltage: 

Electrostatic approach 
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Distribution of charge induced by a gate voltage in a graphene strip is investigated. We calculate 
analytically the charge profile and demonstrate a strong(macroscopic) charge accumulation along 
the boundaries of a micrometers-wide strip. This charge inhomogeneity is especially important in 
the quantum Hall regime where we predict the doubling of the number of edge states and coexistence 
of two different types of such states. Applications to graphene-based nanoelectronics are discussed. 

PACS numbers: 73.63.-b, 81.05.Uw, 73.43.-f 



The new material graphene, a monolayer of carbon 
atoms with honeycomb lattice structure, is attracting a 
lot of interest since 2005 when the first transport mea- 
surements in this material have been reported [H, 0, 9 . 
The interest in 2D electron gases in graphene originates 
from the Dirac-like spectrum of the low-energy quasipar- 
ticles [4] . Several prominent phenomena have been inves- 
tigated in this "relativistic" system both experimentally 
and theoretically, including quantum Hall effect 
[SHI) weak localization and other effects of disorder [_ 
superconducting proximity effects d, Qil HH, 03 , etc. 

In the experiments [J H, 01, mechanically exfoliated 
graphene samples were separated from the metallic gate 
by a b s» 0.3/im wide insulating layer (SzC^)- The 
width of the insulator is dictated by the necessity to iden- 
tify optically the single-layer graphene. In the undoped 
graphene (half filling), the charge of the conduction elec- 
trons is compensated by the charge of the carbon ions 
forming the lattice. By applying a large (V g < 100V) 
voltage V g to the lower gate one induces a considerable 
(rie/Vg ~ 7.2 10 10 cm~ 2 /V [lj) uncompensated charge 
exn e in the graphene plane. This extra charge is screened 
by "image charges" induced in the metallic gate. How- 
ever, since the images are located 0.6fj,m below graphene, 
such a screening becomes effective only in the central re- 
gion of several microns large graphene samples. As a 
result, the charge distribution cannot be homogeneous. 

In this paper, we calculate analytically the charge dis- 
tribution in the graphene strip and demonstrate a strong 
increase of the charge density near the strip edges (nu- 
merically, a charge accumulation near the edges has been 
seen in Ref. [Hj]). For a gate voltage of w 10V the dis- 
tance between the excess electrons in the sheet is of or- 
der ~ lOnm. This means that for the 0.1 4- 1/zm wide 
strips one may speak of a continuous charge distribution 
and determine the latter minimizing the electrostatic en- 
ergy of the electrons. In semiconductor heterostructures 
the electron redistribution has been discussed in the con- 
text of compressible/incompressible QH stripes forma- 
tion [3]. However, in that case electrons were confined 
by a smooth potential, which resulted in a continuous 
charge density profile at the edge. As we will see, at the 
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FIG. 1: Considered experimental setup. Charges (doped 
graphene strip) are placed on the surface (y = 0) of 0.3/im 
thick insulating layer (Si02, dielectric constant e = 3.9) 
above the metallic (n + Si) gate. Image charges (charged 
strips) are shown above and below the graphene plane. 



sharp graphene edge the charge develops a 1 
larity. 
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The charge inhomogeneity discussed in the present pa- 
per develops at the scale ~ 0.1 micrometer. So this is 
a macroscopic effect that should have clear experimental 
consequences. The distribution of classical excess charges 
found below is valid for any metallic strip. However, only 
in graphene the excess charge density coincides with the 
carrier density and determines directly the Fermi mo- 
mentum, pi? a -Jul- The charge accumulation at the 
graphene boundaries is especially important for the quan- 
tum Hall effect where we predict coexistence of two types 
of edge states [ij], [H, EH ■ The strong dependence of the 
charge density on the strip width may have interesting 
nanoelectronic applications, discused at the end of the 
paper. 

Electrostatic potential created by the charge located 
on the surface of the insulator, both above and inside, is 
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FIG. 2: Screening of the electron potential on the surface 
of the insulator (e = 3.9), as described by Eqs. I)2l3[l . Solid 
line shows the screened potential; upper dashed line is the 
potential of the charge on the surface without image charges 
included, Eq. ([1]), and the lower dashed line stands for the 
separate contribution due to image charges. We measure all 
distances in units of the width of insulating layer, b = 0.3/im, 
which is fixed in the experiment. For the in-plane distance 
r = \/x 2 + z 2 > & we have (f> ~ 2b 2 e/e 2 r 3 . 



given by [Vi 
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where the dielectric constant esio 2 = 3-9. In order to de- 
scribe the potential created by the charge placed on a di- 
electric layer (—6 < y < 0) with a metallic gate attached 
underneath (y < —b), one has to consider potentials of 
a string of image charges, as illustrated by Fig. 1. [We 
reserve the coordinates x and z for the graphene plane, 
or the surface of the insulator.] Two different expres- 
sions describe now the potential inside the insulator, at 
-b < y < 0, 
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and above it, at y > 0, 
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Here £ = (1 — e)/(l + e), and the vector b, \b\ = 6 is di- 
rected along the y-axis, perpendicular to graphene plane. 
One may easily show that the potential Eqs. (|2l3p sat- 
isfy the boundary conditions <j)\ — <f>2 and £\d§\jdy — 
£2<902/ dy at the surface of the insulator, y = 0, and 
Eq. ([2]) gives <fi — const = at the metallic surface, 
y = —6. Fig. 2 shows the in-plane potential found from 
these formulas. This potential describes the electron in- 
teraction in mechanically exfoliated graphene. 

Let us consider now a narrow graphene strip with the 
width 2a, such that a <C 6, directed in the x, z plane 
along z-axis. Since the charge is distributed uniformly 
along the strip, the potential both inside and above the 



insulating 5*102 may be obtained from the real part of a 
holomorphic function w(Z) as <$> = Rew(Z), Z = x + iy, 
Aw(Z) = 0. In particular, the function 
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is a solution of the Poisson equation A0o = 
the charge density 
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where 5(y) is the delta- function and a is the charge per 
unit length of the strip. The factor 2/(1 + e) in Eqs. (|1I4|) 
accounts for the polarization of the dielectric substrate. 
The function po ^ is the equilibrium charge distribution, 
since the potential 4>q, is constant on the strip, (j) = at 
— a < x < a,y = §4§. The inverse square root edge 
singularity in Eq. (0, p ~ l/V^ — a, drastically differs 
from the square root density profile p ~ y/x at the soft 
wall edge in the conventional heterostructures [l4[ ■ 

Straightforward generalization of Eq. ([2]) gives the po- 
tential inside the insulating layer sandwiched between the 
metalic gate and narrow (a <C b) graphene strip 

n=o 1 + £ 

Equations ([5]) and ^ are the central result of this pa- 
per describing the charge and potential distributions in 
the narrow mechanically exfoliated graphene strip. Be- 
low we show that these results remain quantitatively ac- 
curate even for sufficiently wide strips, when a « b. 

After the image charges are added the potential on a 
strip acquires a small, ~ (a/6) 2 , coordinate dependent 
correction. From Eq. ^ at y — 0, —a < x < a we find 
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0.175. 



For e = 3.9 we found C » -0.31 and C* 2 ? 

The coordinate dependence of the potential Eqs. (|6l7p 
on the metallic strip should be compensated by a proper 
charge redistribution. Since <fi(x), Eq. ([7]), increases to- 
wards the edges, one should transfer some charges from 
the boundaries to the strip center. To find the equilib- 
rium distribution for finite a/6 we add a series of "mul- 
tipole" corrections to the potential cf>o, Eq. ^ [18| . 
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The same corrections should be added to all image strip 
potentials in Eq. ©. Corresponding corrections to the 
charge density, Eq. ([5]), are 
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FIG. 3: Potential (thick lines) and density (thin lines) across 
the strip of width 2a = 2b — 0.6pm. The dashed lines show 
the curves for the density po, Eq. ([5| and the solid ones - for 
the density p = po — 0.12p2, Eq. ((8j). At the plateau one has 
4>{x) « 0.99a. A gate voltage V g = 100V creates in such a 
strip an averaged electron density (n) = 11.6 x 10 12 cm~ 2 and 
the minimal density n(0) = 8.25 x 10 12 cm~ 2 , while an infinite 
graphene plane gives n<x, = 7.2 x 10 12 cm -2 0]. Semiclassical 
approximation used for the density calculation breaks at the 
distance Sx/b > 0.05(V/V S ) 2/3 from the boundary. 

For example, pi — (2x 2 /a 2 — 1)pq. Still the singularity 
p ~ l/y/x ~ a at the edge is generic for any strip width. 



To compensate the 



term in Eq. ([7]) it is enough 



to consider the first correction only, a-i = — 0.175(a/6) 2 . 
This allows us to approximate the equilibrium distribu- 
tion <fi(x) = const with the accuracy better than 0.2% for 
a < 0.56 (strip width 2a < 0.3pm). A simple formula 
(both V g and a have dimensionality charge/ distance) 



V g = cr[0.82m(6/a) + 0.88 + 0.29(a/6) 2 



(10) 



relates in this case the linear charge density a in the 
narrow strip to the applied gate voltage. 

An appropriate fit with only two parameters «2 , «4 7^ 
allows us to reach <f)(x) w const on the strip with the 
accuracy ~ 0.5% even for a = 56 (strip width 2a — 3/xm). 
Remarkably, even for such a wide strip the amplitude of 
the singularity is reduced only by a factor 0.55 

compared to the simple formula, Eq. |(SJ). Fig. 3 shows 
results of the single parameter fit, 0,2 / 0, for a = b. 

The classical charges equilibration condition, applica- 
ble for any metallic strip, allows to find the inhomoge- 
neous electron density, Eqs. (|5I9[) . but leads to the con- 
stant potential in plane. A nontrivial graphene-specific 
potential profile across the strip appears due to quantum 
effects. Quantum dynamics of electrons in graphene is 
described by the Dirac equation 



[vf(t x Px ± T y p.y) +U(x)]il>± = sip±, 



(11) 



where r x ,y are the Pauli matrices interchanging the sub- 
lattice index on the honeycomb lattice. [Strictly speaking 
we should write here (r x p x ± r y p z ), since we use coordi- 
nates x,z for the graphene plane, not x,y used usually 
in the literature.] The two signs ± correspond to two 
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FIG. 4: Schematic charge distribution p(x) in units of 
10 11 cm -2 in the QH regime for a strip width 2a — 2b = 
0.6/xm, B = 10T, V g = 5V (thick solid line). The ~ Ij^fx 
in crease of density at the edge is stopped at Aa; ~ Ib = 
\JhcjeB ~ 8nm. Short-dashed line shows the electrostatic 
solution, the same as in Fig. 2. We assume valley-degenerate 
Landau levels, and choose the Zeeman splitting -Ezecman = 
0.25i5o [6J. Lower dashed curves show the effective potential 
for different Landau levels U e g(x) = U (x) + y~N Eo ± -Ezeeman 
(in units of 0.2eV). All electronic states with U e g < 0(> 0) 
are occupied(empty). Regions with U e s = correspond to 
partially occupied Landau levels (compressible stripes). The 
figure shows coexistence of two types of edge states: Com- 
pressible stripes in the center and usual noninteracting edge 
states at the borders. 

valleys in graphene and ip± are envelope functions. So- 
lutions of Eq. (fTTj) are double degenerate due to spin 
and vf ~ 10 8 cto/s. The Pauli principle prevents all un- 
compensated electrons in the strip from having the same 
zero momentum, p — 0, as was assumed in the electro- 
static solution, Eq. (O. To account for the coordinate 
dependent electron density we introduced in Eq. (|11[) a 
potential U(x) < 0, while keeping the zero Fermi en- 
ergy, Ep = 0. For a large charge density, the potential 
U(x) varies slowly on the scale of the wave length A, 
which allows us to introduce the local Fermi momentum 
Pf = H/Xp = \U(x)\/vf- The density of electrons can be 
found in the Thomas-Fermi approximation as 
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The 2-dimensional density of electrons n e is related to 
the 3-dimensional charge density used in Eqs. (|5l9p as 
p = en e 5(y). Thus for narrow strip, a <C 6, we find from 

Eqs. (isira 



U(x) 



x 2 ) 



-1/4 



(13) 



For the gate voltage V g — 100V and the strip width 2a = 
0.6pm we estimate (7(0) = — 0.335eV. This quantum 
(U(x) (x h) correction to the electrostatic potential on the 
strip describes locally the position of the Dirac crossing 
point with respect to Fermi energy. 

The semiclassical approximation used here is justified 
provided \dXp/dx\ <C 1. Thus we may use Eqs. (J3J 




FIG. 5: Five-terminal device for measuring separately the 
Hall currents in inner and outer edge states. Only the current 
carrying channels are shown. The two currents cancel each 
other in the gate 2, while the split gates 4 and 5 measure each 
the current due to several inner or outer states. 

only at distances Sx > a 1 ^ 3 (e/cr) 2 / 3 from the strip edge. 
At Sx ~ a 1 / 3 (e/cr) 2 / 3 the singular increase of both the 
density Eq. © and the potential Eq. IT51) is stopped in 
a way dependent on the details of the graphene edge. In 
particular, the maximal value of the density is n max = 
const x (cr 2 /e 2 a) 2 / 3 with const ~ 1 depending on the type 
of the edge. 

Experimentally, the conductivity £ of graphene in- 
creases linearly with the gate voltage [H, 12, I3- This 
implies £ cx vfPf oc ^fp. The increase of the carrier 
density near the edges of the strip should lead to an in- 
homogeneous current density distribution j (x) oc W p(x). 
Presence of a (moderately strong) disorder should not 
change the density distribution in the strip. 

The non-monotonic charge distribution across the 
graphene strip (l9l | should be especially important in the 
QH regime when the electron transport is due to exis- 
tence of edge states [15j. The charge density in this case 
is roughly given again by Eq. ([5]) with the 1 / \fx edge sin- 
gularity smoothed at the distances ~ Ib = \JhcjeB. The 
number of occupied Landau levels as a function of the 
transverse coordinate x first increases almost abruptly 
(at the length ~ l B ) at the strip edge and then decreases 
towards the strip center (length scale ~ a). This leads 
to formation of a double set of QH edge states, having 
a very different microscopic nature (see Fig. 4). At the 
graphene edge one may effectively neglect the electron 
interaction and consider the QH edge states formed by 
the last occupied electron state on the n-th branch of so- 
lutions of single-particle Dirac equation E n (xo) (the n-th 
Landau level) [Hill- 

This is an adequate description 
of the outer edge states in graphene @ . Away from the 
boundaries the electron repulsion transforms the narrow 
(~ Ib) edge states into the wide compressible stripes with 
a macroscopic number of electrons having similar ener- 
gies U e s(x) = const. These compressible stripes alternate 
with the incompressible ones having a constant electron 
density n e (x) — const [T3 |. 

The energy of the N-th Landau level in graphene is [5| 
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FIG. 6: Creation of quantum dots via the charge accumula- 
tion in narrow constrictions, a). Single QD. b). Double (par- 
allel) QDs. c). Side coupled QD. Shaded areas show the lakes 
of large electron density. A constriction in biased graphene 
strip works not as a Quantum Point Contact, but as a QD. 

component. So, we find the number of occupied Landau 
levels at a given point in a strip 



N(x) = (U(x)/E ) 2 = n e hc/4eB. 



(15) 



E w VnE , E = hv FX /2eB/hc. 



(14) 



There is a series of such levels for each spin and valley 



This formula shows the smoothed number of occu- 
pied Landau levels for N(x) ^> 1. Going beyond 
this approximation reveals the compressible and incom- 
pressible striped QH states shown in Fig. 4. The 
picture is schematic since the details of compress- 
ible/incompressible stripes are not described by the 
smooth Eq. (fT5"l) . [Still this may be done electrostati- 
cally, see Ref. [14 |, in the regions there dU jdx <C Eq/Ib-] 
The form of the density close to the edge, Sx ~ Ib, as 
well as the physics of outer edge states ja, Hcj , depends 
on the form of graphene edge. Nevertheless we may say 
that the lower gate voltage V g « 5V should be sufficient 
to create several edge states of both kinds for the strip 
width 2a « 0.6pm. 

Each QH edge channel supports the electric current 
flowing only in one direction. For the appropriate sign of 
the bias voltage V , the channel may carry the current j = 
e 2 V/h. Since electrons belonging to inner and outer QH 
edge channels drift in opposite directions, corresponding 
currents have a tendency to compensate each other. In 
Fig. 5 we suggest a simple five-terminal device, which 
would allow to measure separately the currents carried by 
the outer and inner channels [2lj . The QH effect in such 
a setup would be seen at parametrically smaller values of 
the gate voltage than in the existing experiments 0, Q ■ 

A striking consequence of the result, Eq. (J7J), (see also 
Ref. [l3|]) for the strip potential is that for a <C b the in- 
plane charge density is inversely proportional to the strip 
width. (The linear density of charge a — J p(x)dx de- 
pends only logarithmically on the strip width, (jTUJ) , hence 
p(x) ~ 1/a). This offers a possibility of creating lakes of 
a large charge density, quantum dots (QD), by cutting 
narrow constrictions in the graphene strip. Examples of 
such devices are shown in Fig. 6. This semi-mechanical 
way of confining electrons (potentials that appear due 
to the strip narrowing lead to the longitudinal confine- 
ment) may be complementary to the pure electrical way 
of fabrication of QDs in graphene 0, HH, IH| • 

Without the electric field doping (V g = 0) graphene 
behaves as a hole metal [l|. The shift of the Fermi en- 
ergy away from the Dirac crossing point is attributed to 
an unintentional doping of the film by absorbed water. It 
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is compensated by application of a sufficient lower gate 
voltage (~ 40eV [lj). However, as we have shown, the 
gate-induced charging is nonuniform, and it is impossible 
by varying V g to reach the Dirac point simultaneously in 
the whole sample. This charge distribution evolves dif- 
ferently for V g below and above the value corresponding 
to the strip minimum conductivity (for example, the lo- 
cal crossing of the Fermi energy by the Dirac point is 
shifted towards the center or the boundary of the strip) . 
This may explain the I-V characteristics assymetry in 
graphene (observed e.g. in Ref. [l2|). 

In conclusion, in this paper we predict and describe 
the macroscopic charge accumulation along the bound- 
aries of graphene strips, made of experimentally used me- 
chanically exfoliated films, for moderate (< 10V) lower 



gate voltages. Information about the local Fermi mo- 
mentum and charge density pf oc y/rQ may be extracted 
from the STM measurements of the density of states in 
graphene 25|. The average charge density (p) = a /2a 
for given gate voltage also strongly increases for narrow 
strips (< 0.5/im) as described by Eq. (|10|) . Transport in 
graphene would be especially sensitive to the predicted 
charge accumulation in the experiments 0, 3 in the QH 
regime, where the two kinds of edge states [Til [HI, [Hj] 
should coexist in the same sample. Experimental setup 
capable to measure currents carried by different edge 
states is suggested (Fig. 5). 
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